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ABSTRACT 



Frequencies of magnetic patch processes on supergranule boundary, namely 
flux emergence, splitting, merging, and cancellation, are investigated through an 
automatic detection. We use a set of line of sight magnetograms taken by the 
Solar Optical Telescope (SOT) on board Hinode satellite. We found 1636 positive 
patches and 1637 negative patches in the data set, whose time duration is 3.5 
hours and field of view is 112" x 112". Total numbers of magnetic processes are 
followed: 493 positive and 482 negative splittings, 536 positive and 535 negative 
mergings, 86 cancellations, and 3 emergences. Total numbers of emergence and 
cancellation are significantly smaller than those of splitting and merging. Further, 
frequency dependences of merging and splitting processes on flux content are 
investigated. Merging has a weak dependence on flux content only with a power- 
law index of 0.28. Timescale for splitting is found to be independent of parent flux 
content before splitting, which corresponds to ~33 minutes. It is also found that 
patches split into any flux contents with a same probability. This splitting has a 
power-law distribution of flux content with an index of —2 as a time independent 
solution. These results support that the frequency distribution of flux content in 
the analyzed flux range is rapidly maintained by merging and splitting, namely 
surface processes. We suggest a model for frequency distributions of cancellation 
and emergence based on this idea. 

Subject headings: Sun: granulation — Sun: photosphere — Sun: surface magnetism 
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1. Introduction 

How magnetic structure on the solar surface is constructed and maintained is one of 
the fundamental issues in solar magnetic field observation. It is important for the statistical 
understanding of solar activities because they are triggered by magnetic activities on the 
solar surface. It may also give a quantitative restriction to solar dynamo problem. One 
approach is to investigate a frequency distribution of magnetic flux conte nt on the solar 



surface. Some author s found an exponential distribution of flux content ( iSchrijver et al. 



1997 




Haeenaar 


1999 


). On the o 


;her hand, other authors 


(Wane et al. 


1995; 


Parnell et al. 


2009; 


Zhang et al. 


2010 


)• 



bund a power - law d istribution 



Parnell et al 



( 1200a ) reported a 



power-law distribution with an index of —1.85 between 2 x 10 17 Mx and 10 23 Mx, which 
means that magnetic patches from large active regions to small patches in quiet network 
are described by a single flux distribution. They suggest the idea that either all surface 
magnetic features are generated by the same mechanism or that they are dominated by the 
surface processes. 

The next arising question is how the frequency distribution is achieved and sustained. 
Magnetic processes, namely flux emergence, splitting, merging, and cancellation of magnetic 
patches on the photosphere are thought to change and maintain the frequency distribution. 
Relationship between these mag netic processes and t he flux distribution is described 



by magneto-chemistry equation (iSchrijver et al 



19971 ). Based on this equation and 



some assumptions, they found a time ind ependent solu tion for an exponential frequency 



distribution of flux content. Furthermore 



Parnelll ( 120021 ) found a particular solution for an 



arbitrary frequency distribution of flux content. Both solutions assume re-appearances of 
canceled fluxes and a detailed balance between any two fluxes, which should be verified 
observationally. 



From the view of flux balance on the photosphere, flux emergence and cancellation are 
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especially investigated because they have a direct relation with a flux exchange through the 
photosphere. Flux emergence, which is observed as a divergence of opposite polarities in 
magnetogram, is thought to be a flux ascension from below the photosphere. It produces 
flux increases of both polarities in line-of-sight mag netograms . The frequency distribution 



of emerging flux is investigated by several authors. 



Hagenaarl (120011 ) fo und an exponential 



distribution by using full disk magnetograms obtained by SOHO/MDI. Thornton & Parnell 



( 120 111 ) found a power-law distribution from active region to inter-network field by 



using Hinode/SOT magnetogram. Flux cancellation is defined as a convergence and 



Martin et al. 


1985; 


Livi et al. 


1985) 



fluxes on the s olar surface. Two physical models, U-loop emergence and Q 



are proposed (IZwaan 



Yurchvshyn &: Wang 



2001 



1987) and many autho r s tried to distingu i sh them 



Chae et al 



2002 



loop submergence 



Kubo et al.ll2010t llida et al. 



arvev et al. 



2010 



1999; 



Chae et al 



20101 ). Flux replacement timescales by emergence and cancellation are physical quantities 



representing importan ce of these processes on the solar sur 



days to several hours (ILivi et al. 



1985 



Schrijver et al. 



1998 



ace. But it varies from several 



Hagenaar 



20011 ) and it is not 



clear why this timescale decreases as the spatial resolution becomes higher. 

Merging and splitting should also take important roles in flux maintenance because 
magnetic patches change their flux content through these processes. The action of convective 
motion should differ among patches with different flux content. However, there are fewer 
reports for merging and splitting compared to those of emergence and cancellation. The 
role of these processes in flux balance still remains unclear. 

We investigate the frequencies of magnetic patch processes, namely flux emergence, 
splitting, merging, and cancellation based on observations in this series of papers. This first 
paper is mainly devoted to investigations of total flux amount of magnetic processes and 
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frequency distributions on flux content of merging and splitting. The purpose of this paper 
is to clarify what the processes dominate the frequency distribution of flux content. Our 
final goal of this series of papers is to understand how the flux distribution is maintained on 
the solar surface. We concent r ate magnetic process es on the network because most flux is 



contained there (Martin 



1990 



Schrijver et al. 



19981 ). The detail of data set is explained in 



section 2. We explain our definition of a magnetic patch and magnetic processes in section 
3. The results are shown in section 4 and related discussions in section 5. 



Observation 



We use line-of-sight magnetograms near the disk center obtained by N arrowband Fi 



Imager (NFI) of So 



2007 



Tsuneta et al- 



ar Optical Telescope 



2008 



Ichimoto et al. 



(SO 



2008 



onboard Hinode satellite (IKosugi et al 



Shimizu et al 



2008 



Suematsu et al 



ter 



2008|) 



We need to use Hinode/ NFI magnetograms in this study bec ause the detection 



Parnell et al. 



imit of 



(120091 )) and we 



MDI high-resolution data is ~ 10 18 Mx (see a Figure 5 (a) by 
can not set a flux range wide enough for an analysis of the network field. Time period of 
the data set is from 0:33UT to 4:08UT on 2009 November 11. NFI observed magnetograms 
near disk center by using Na I Dx resonance line at 5896 A during this period. Full field of 
view is 112" x 112". Figure 1 shows an example of magnetograms used in this study after 
the pre-process explained in the next section. Some network cells, which has a typical size 
of 20 — 40", are included in the field of view. Total number of magnetograms is 199 during 
the whole observational period. Time interval between each magnetogram is 1 minute. The 
interval is sometimes 2 minutes due to the lack of data. The region is a quiet Sun because 
there are no active regions on the solar disk during the observational period. 

Some corrections of SOT magnetogram are done before detecting magnetic processes. 
We use f g_prep.pro procedure in SolarSoftWare (SSW) package for a correction of dark 
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current and flat field of CCD camera. The data is rotated to the position at November 11 
2:03UT when the region is near the disk center by using drot_map. pro procedure in SSW 



2010h by 



package. We remove a columnwise median offset of CCD camera (lLamb et al. 
subtracting a median value of each column from all pixels in the same column. Small spatial 
fluctuation of field of view is removed by making a correlation between consecutive images. 
Then the magnetograms are averaged over 3 continuous images and 3 pixels for smoothing. 

The observed circular polarization signals (CP) are converted to the line-of-site 
magnetic field strength by the following method. We multiply CP by a conversion coefficient 
to obtain actual magnetic field strength in the next step. This conversion coefficient 
is determined by comparing the SoHO/MDI HR magnetogram and CP image on 2009 
November 11 1:39UT. We spatially smear the SOT data to two times of MDI pixel size 
(R. A. Shine in private communication) and make a linear fitting between CP in SOT and 
magnetic field in MDI. Because there are too weak or strong pixels in SOT data which are 
out of linear range, we make the fitting in the range from 30 G to 100 G as a magnetic field 
strength. Figure 2 shows a scatter plot of SOT data number and magnetic field strength 
obtained by MDI. The solid line shows a result of linear fitting. We obtain 9067.98 G/DN 
as a conversion coefficient. 



3. Identification of Magnetic Processes 

We define magnetic processes and explain our method for the identification of them in 
this section. Figure 3 represents schematic pictures of four magnetic processes detected in 
this study and Figure 4 summarizes our definition of magnetic processes. Note that our 
definitions of these processes are valid even when more than two patches are involved in the 
process but invalid when processes which increase and decrease flux content occur at the 
same time. 
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3.1. Detection and Tracking of Magnetic Patches 



We use a clumping method for a detection of magnetic patches: each patch is picked 



up as a clump of marked pixels having magnetic strength beyond a given threshold ( IParnell 



20021 ). The adopted threshold is one sigma obtained by fitting the histogram of the signed 
magnetic stren gth in each magnetog ram, i.e. o ~ 5 G pixel ~ 6.8 x 10 14 Mx. This value is 



close to that of 



Parnell et al. 



(120091 ) . A magnetic patch in a clumping method corresponds 



to a massif of magnetic patches in a downhill method an d a curvature method which are 



used in some previous studies (T lWelsch fc Longcope 



20031 ; ?). We pick up magnetic patches 



with sizes beyond 81 pixels for focusing our analysis on the network magnetic field. The 
validity of this choice is demonstrated in Figure 5, in which the network structure is seen 
when adopting the 81-pixels threshold. 

We track the motion of magnetic patches in consecutive images after the detection of 
them. Patches are mar ked as identical when they have a spatial overlapping in continuous 



images (jHagenaar 



19991 ). It should be noted that a travel distance of magnetic patches in the 
data interval (~lmin) is nearly up to one pixel size, namely 2 km s _1 x 1 min = 120 km ~ 1 
pixel. In high-resolution magnetogram, more than one patch in a previous image often 
have spatial overlappings with one patch in a consecutive image and vice versa. To clear 
up this problem, we set two conditions when tracking patches. First we investigate spatial 
overlappings of patches from those with larger flux contents. This is based on a concept 
that a smaller patch has a greater tendency to fall below a detection limit of the analysis by 
splitting and cancellation. Second we select a patch with the most proximate flux content 
in case of overlappings of more than one patch. Tracking paths of detected patches become 
unique with these conditions. 
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3.2. Merging and Splitting 

Merging is a process where more than one patch of same polarity converge and 
coalesce to one patch (Figure 3a). We define a merging as an event with two conditions. 
A merging event is defined as a feature which satisfies the conditions 1) that there are 
two or more parent patches in a previous magnetogram overlapping one daughter patch 
in the consecutive magnetogram, 2) and that more than one of the concerned patches in 
the previous magnetogram disappear in the time interval. Panels in the first row of Figure 
6 show an example of detected mergings. Two patches converge in the first and second 
images. They coalesce to one patch in the third image. 

Splitting is a process where a single patch is divided into more than one patch (Figure 
3b). A splitting event is defined as a feature which satisfies the two conditions 1) that there 
are one or more daughter patches in a magnetogram overlapping one parent patch in the 
previous magnetogram, 2) and that more than one of the concerned patches in the latter 
magnetogram appear in the time interval. Panels in the second row of Figure 6 show an 
example of detected splittings. One patch stays in the first and the second panel. It splits 
into two patches between the second and third panel. 

3.3. Emergence and Cancellation 

The detailed explanation for the detection technique of emergence and cancellation 
will be given in the next paper in which much more events are detected by using another 
data set, making it possible for us to conduct a statistical study of them. Emergence and 
cancellation are defined as a pair of flux change events of each patch in opposite polarities 
within a certain distance. We pick up not only complete cancellations but also partial 
cancellations in this definition. 
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Panels in the third row of Figure 6 show an example of emergences. In the second 
panel, positive and negative patches appear but the positive patch is already coalesced with 
a pre-existing patch. These emerging patches continue to separate against each other and 
become distant in the last panel. Panels in the last row of Figure 6 show an example of 
partial cancellations. Opposite polarities are converging in the first and second panel. They 
start canceling between the second and third panel. The negative patch totally disappears 
in the last panel but positive patch remains after the cancellation in this example. 

4. Results 

4.1. General Description 

Along with the total numbers of patches and their flux density averaged over field of 
view and the observational period (in Mx cm" 2 ) and the frequencies of events in magnetic 
flux density (in Mx cm" 2 s" 1 ) are summarized in Table 1. Total numbers of detected 
negative and positive patches are 1636 and 1637, respectively, enough for a statistical study. 
As for magnetic processes, merging and splitting are much more frequent than emergence 
and cancellation. Total numbers of merging for positive and negative patches are 493 and 
482. Total numbers of splitting for positive and negative patches are 536 and 535. These 
are enough for a statistical study. The number of emergence is 3 and that of cancellation is 
86, which are not enough for a statistical study. 

To determine what magnetic processes dominate a flux balance in this region, we 
investigate total flux amounts of patches and those related to the magnetic processes per 
unit area per unit time. The results are shown in Table 1. The averaged flux density are 
2.53 Mx cm -2 and 3.60 Mx cm -2 for positive and negative patches respectively. They 
correspond to the total flux amount of 1.63 x 10 20 Mx and 2.32 x 10 20 Mx for our field of 
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view. The rate of averaged flux density involved in merging processes is 1.65 (3.53) x 10~ 3 
Mx cm -2 s _1 for positive (negative) patches. We define flux amount of merging process 
as that of parent patches. Flux replacement timescale by positive (negative) merging, 
T mrg = $tot/(<9$tot/<9£| mrg ), is evaluated as 1.53 (1.02) x 10 3 sec from these values, which 
is much shorter than the estimat ed replacement t i mescale by cancellation and emergen ce 



reported in the previous studies (ILivi et al. 



1985 



Thornton & Parnell 



Schrijver et al. 



1998 



Hagenaar 



2001 



20111 ). The rate of averaged flux density involved in splitting processes 



is 1.48 (3.03) x 10 -3 Mx cm -2 s _1 for positive (negative) patches. We define flux amount 
of splitting process as sum of those of parent patches. In the same manner as merging, flux 
replacement timescale by positive (negative) splitting is evaluated as 1.71 (1.19) x 10 3 sec. 



4.2. Frequency distribution of flux content 



Figure 7 shows a frequency distribution of flux content. Red and blue dashed histogram 
denote those of positive and negative polarities, respectively. Solid histogram denotes sum 
of them. Black dashed line indicates the fitting result of it. The distribution is dropping 
down below 0th = 10 17 ' 5 Mx, which is suggested as a detection limit in this study, below 
which only a limited number of the patches are detected. Dashed line is a fitted power-law 
result between 10 17,5 Mx and 10 19 Mx. We make a least-square fitting of power-law function 
with a form, n{4>) = n (<p/<j) )-~> . We obtain n = 2.42 x 10~ 36 Mx" 1 cm" 2 , O = 1.0 x 10 18 
Mx, and 7 = 1.78 ± 0.05. The error me ans one sigma error. The fitted power-law index is 



sensitive to fitting range as reported by 



Parnell et al. 



(120091 ) . The index varies from —1.78 



to —1.91 with a minimum fitting range of 10 17 5 to 10 18 Mx. Total error of our fitting of 
the power-law index bec omes 0.18 with considering this error. This result is consistent 



with 



Parnell et al 



(120091 ) . which reports —1.85 ± 0.14 as an power-law index of a frequency 



distribution of flux content. 
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4.3. Probability distributions on flux content of merging and splitting 

We investigate frequency distributions on flux content of merging and splitting for one 
patch, which are defined as frequency distributions of processes divided by a frequency 
distribution of flux content. The value of these distributions represent probability of the 
processes. Further, the detection limit t h is taken into the account in this study. We call 
them apparent probability disributions of processes. 

Figure 8 shows the apparent probability distribution of merging. We make a 
least-square fitting in a range of 10 17 5 - 10 19 Mx, where number of merging event is enough 
for fitting. The fitting form is 



where po.mrg is a reference probability, 0o is a reference flux content, and /3 mrg is a power-law 
index of probability distribution of merging. We obtained po.mrg = (2.52 ± 0.08) x 10~ 4 
s -1 and /3 mrg = 0.28 ± 0.05 for positive patches, p ,mr g = (2.52. ± 0.08) x 10~ 4 s -1 and 
Anrg = 0.26 ± 0.05 for negative patches, and po,mr g = (5.12 ±0.11) x 10~ 4 s -1 and 
Anrg — 0.28 ± 0.04 for both patches with O — 10 18 Mx. The errors mean one sigma error of 
least-square fitting. 

Figure 9 shows the apparent probability distribution of splitting. The strong increase 
in the range larger than 10 19 Mx is caused by lack of the patch number in the analysis. 
On the other hand, there is a drop in a flux range near th where the number of patches 
is enough for a statistical study. We interpret this dropping as an effect of splitting into 
the area below t h- This effect is evaluated in the discussion in Section 5.1. We see that 
probability of splitting is almost constant as 1.0 x 10~ 3 sec -1 , which means a timescale 
of 33 minutes, in the range enough above t h, 3.0 x 10 18 — 1.0 x 10 19 Mx. It means that 
frequency of splitting is independent of parents' flux content. 




(1) 
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5. Discussion 

5.1. Probability density distributions of merging and splitting 

We evaluate probability density distributions of merging and splitting terms in magneto- 
chemistry equation from the observational results. Because probability distributions, which 
we obtained in this study are obtained by integrating them on flux content once, we have 
to put at least one assumption to evaluate them. 

The merging function l(x,y) is given as follows. We obtain the form as a probability 
distribution from the definition of l(x,y) as 



npAPP 
u mrg 

dt 



POO 

/ n(x)l((t), x) d(j). (2) 

J 4>th 



'<f>th 

By comparing with the observational result, we obtain 

roo / ^ x /3 mrg 

/ n(x)l((j),x)dx =J9 ,mrg [ ~T ) (0 > 0th)- (3) 

J<t>th \0O/ 

In the left-hand-side of this equation, the variable <fi appears only in the l(<p,x). We assume 
a simple form satisfying this relation namely, 

l(x,y) oc x 13 ^. (4) 

From the symmetry of l(x,y) = l(y,x), this relation deduces 

Substituting it and observational result of n((f>) into Eq.(3), we obtain 

k = ^ • (6) 

^o"^ 6 / X^^dx 
J<pth 

The upper value of the integration is limited on the actual Sun and we put it as max - We 
substitute the value obtained in our study, namely n = 1.21 x 10~ 36 Mx" 1 cm -2 , 7 = 1.78, 
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Po,mr g = 2.56 x 10- 4 s-\ /3 mrg = 0.28, max = 10 19 Mx, th = 10 17 - 5 Mx, and O = 10 18 Mx, 
and obtain l = 7.24 x 10 13 cm 2 s -1 . 

The splitting function k(x, y) is given as follows. From our observations, the probability 
distribution of splitting events dP^f^ ((f)) /dt is suggested to be independent of the parent 
patch flux: 

^pAPP 

— t^— (0) = k = constant (for all 0). (7) 

C6 

This claim is observationally supported at least in the range > t h (Figure 9). The 
drop off below th is discussed immediately below. If the splitting ratio between daughter 
patches is randomly determined, i.e. 

d 

— [k(x,(f>-x)] = (for < x < 0), (8) 



then we obtain, 



k(x,y) = -^. (9) 
x + y 



When the flux content of the daughter patch is below t h, such events are not recognized as 
a splitting event in our procedure. The probability distribution will be given as 



npAPP r<t>-<f>th o fk , 

-%r-W = / 4>-x)d* = fc o(i - =^). (io) 

dt J fa 



Black and blue dashed curves in Figure 9 indicate analytical curves with fc = 1.0 x 10 -3 
s _1 and k = 5.0 x 10~ 4 s -1 respectively. This curve fits the drop of the observational line 
well, which supports the above assumptions. We obtain k(x,y) = k /(x + y) as splitting 
function, where ko = 5.0 x 10~ 4 s _1 . 



5.2. The time-independent solution of splitting process 



Since our observations show that merging and splitting are much more frequent than 
emergence and cancellation, it suggests that the former two have much influence on the 
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maintenance on the power-law distribution. We show the time-independent solution by 
splitting although we have not found the time-independent solution by merging and splitting 
The magneto-chemistry equation is consist of source (emergence) term, m erging terms, 



splitti ng terms, and cancellation terms (see the right-hand-side of Eq.(3) in 



Schrijver et al. 



( 119971 )). The frequency of emergence, merge, splitting, and cancellation are represented 
by S((p), l(x,y), k(x,y), and m(x,y) respectively. We use the magneto-chemistry equation 
only including the splitting terms by setting S((fi) = 0, l(x,y) = 0, m(x,y) = 0, namely 

dn(4>) 



dt 



2 / n(x)k(4>, x — 4>) dx — / n((f))k(x, <fi — x) dx 
><!> Jo 



(11) 



where n(4>) is a frequency distribution of flux content 



After substituting k(x,y) = ko/(x + y) into Eq.(ll) and differentiating with 0, we 
obtain 

d 2 n(<f)) k d 



n 



12) 



d<j)dt cj) 2 dfy 

This equation has a time-independent solution n{<ft) oc <p~ 2 . This power-law index of flux 
content is in good agreement with the observational result. This scale free distribution comes 
from constancies of splitting, namely that splitting has a constant timescale independent of 
flux content and constant probability of splitting to flux content. These constancies may 
come from convection dominating patch stability or flux tube instability with a constant 
timescale. We needs further theoretical and observational studies to determine which of the 
hypotheses is the actual scenario on the solar surface. 



1 Terms in t he rig ht-handside of this equation is different from those of Eq.(3) in 



Schrijver et al. 



( 119971 ) at the point of th e 



paper and n(x)k(<f>,x) in 



Schrijver et al. 



split ting from 0, n amely n(x)k(4>, x — 4>) in this 



3). 



(19971). The term in 



Schrijver et al. 



(119971 ) should 



be a typo because we have to multiply the pre-splitting number density here. 
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5.3. Relationship among frequency distribution of flux content, cancellation, 

and emergence 

We suggest a model of relationship among frequency distribution of flux content, 
cancellation, and emergence. The important hypothesis in this model is that the frequency 
distribution of flux content is rapidly maintained regardless of cancellation and emergence, 
which is supported by this study. Figure 10 shows a schematic view of this model. We put 
a power-law distribution of flux content in unit of patches Mx" 1 cm -2 as 



n((/>) = n o[ — 



(13) 



where 0o is a reference value of flux content and no is a reference frequency of flux content. 
The po wer-law index, 7, is derived as 1.5 < 7 < 2 by our observation and the previous 



studies (IParnell et al 



2009 



Zhang et al. 



20101 ). The maximum flux content in the system 
(0max) is assumed to be much larger than the minimum ((/> min ) in the following discussion. 
We calculate N(4>), a total patch number with flux content larger than 0, by integrating a 
flux distribution from <fi to max as 

-7+1 



N{<f>) 



n(4>') dip' 



n o o / 



7-1 \<f>, 



(14) 



Schrijver et al. 



( 119971 ) evaluated a collision rate of opposite patches from a total patch 
number density with assumptions of a constant velocity and a randomness of patch motion 
along network. They obtained the collision frequency, u, as 



_vq_ 



■N 



(15) 



where vq, p, and N t mean a typical velocity of patches, a number density of network cell, 
and a number density of patches respectively. We multiply 1/2 taking the double counting 
into the account. We obtained that the frequencies of merging and splitting are larger than 
that of cancellation in this study. It can be deduced that the frequency distribution is 
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maintained rapidly by merging and splitting compared to the timescale of cancellation. This 
enables us to treat the number density of patches is time-independent in the evaluation of 
cancellation and apply the same analogy to the number density expanded in the dimension 
of flux content. The collision frequency ,dN(<j))/dt\ col , is evaluated as 

dN{<t>) 2 ^ 2 /^^+ 2 



dt 



col 



^0^000 ( _0 

'4( 7 -1)Vp Wo 



(16) 



Note that this total collision number becomes time-independent with an assumption 
of maintenance of a power-law flux distribution. We assume that a total number of 
cancellation, dN(<fi)/dt\ cnc , equals to a total number of collision events of opposite 
polarities, <9iV(0)/<9t| col , namely 



dN(cf>) 



dt 



dN(<P) 



dt 



(17) 



col 



The frequency distribution of cancellation, dn(4>) / dt\ cnc is given by differentiating Eq.(16) 
with 0, namely 



dn(4>) 



dt 



d 2 N((f>) 



d(j)dt 



-27+I 



(18) 



2( 7 -1)VpV0o, 

This assumption means that there is no patches passing through the patches of opposite 
polarity once they collide, which is difficult to check and we justify this assumption from 
the comparison of obtained frequency distribution of emergence in this model and that in 
the observation. We also evaluate a frequency distribution of emergence with assumptions 
of small amount of flux supply from the outside of the system and re-emergences of 
canceled fluxes. These assumptions lead to the relationship that a frequency distribution of 
emergence nearly equates that of cancellation, 



dn(4>) 



dt 



dn(4>) 



emrg 



dt 



v n 2 Q <i)o 



-27+I 



ic 2( 7 -l)v/pV0o 

We compare the power-law index and absolute value of frequency distribution of 
emergence in this model with the observational results. Based on the above discussion 



(19) 
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our observation suggest a frequency distribution of emergence as dn/dt\ emrg ~ 
4. 3 x 10~ 35 x (0/1.0 x 10 16 Mx)" 2 6 Mx" 1 c m~ 2 s" 1 where we adopted p = 1.0 x 10~ 19 cm" 2 



( IHagenaar et al. 



19971 ) and w = 2.0 km s _1 . 



Thornton Sz Parnelll (1201 if ) reports a power-law 
frequency of emergence as dn/dt\ cmrg ~ 2.5 x 10~ 35 x (0/1.0 x 10 16 Mx)~ 2 ' 7 Mx -1 cm -2 
s _1 . The steepness of the power- law distribution and the absolute value in our model are in 
good agreement with the observational results. 

We calculate flux replacement time by cancellation and emergence. The total flux 
amount in the system is calculated from Eq.(13) as 



tot 



(20) 



Since 7 < 2, this result shows that the total flux is dominated by patches with larger 
flux content. On the other hand, we obtain a total flux loss amount by cancellation, 
<9$tot(0)/^| cnc , from Eq.(18) as 



<9$ 



tot 



dt 



cnc •> 0mi 



, dn(cj)') 



dt 



-27+3 



2(7-1X27-3)^ v 

The total flux supply by recycled emergence is also evaluated in the same manner from 
Eq.(19) as 



(21) 



tot 



dt 



emrg •> <t> 



Ll dn{<l>') 



dt 



2^3 




emrg 



2( 7 -l)(2 7 -3)Vp V 0o 



-27+3 



(22) 



Then flux replacement time is evaluated as 



Replace = $tot/( ^tot/^| cnc ) = $tot/( ^tot/^| 

2( 7 -l)(2 7 -3) /0 m£ ^ 7+2 



emrg/ 
27-3 



(23) 



(2 - 7)^0^000^ V 0o / V v>o 
This result is qualitatively consistent wi th the previous result that flux replacement time 



becom es shorter with higher resolution ( IMartin et al. 
200lk 



1985 



Schrijver et al. 



1998 



Hagenaar 
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We summarize our interpretation from the discussion. In our interpretation, a 
power-law frequency distribution of flux content is rapidly maintained by merging and 
splitting, which is supported by the result of a comparison of change rate of flux amount 
by each process. Addition to this, we have found that a power-law frequency distribution 
with an index of —2 is a time-independent solution of splitting. Cancellation is caused by 
convective dominant motion and frequency distribution of cancellation should naturally 
become a steep power-law distribution. Emergence is interpreted as re-emergence of 
submerged flu x by cancellation, which i s consistent with a steep frequency distribution 



of emergence (IThornton &: Parnell 



2011 



j). This should be one of the possibilities but the 
important point of this model is that the apparent flux transport through the photospheric 
layer becomes much more drastic when investigating emergence and cancellati on although 



i njecte d flux amount from the deeper layer is small or even zero. Recent paper, 



Meyer et al. 



(120 111 ), reports their numerical simulations of the magnetic carpet on the photosphere. 
They show that a power-law distribution of flux content is maintained with an input of a 
steep power-law frequency of emergence. The induced frequency of cancellation becomes a 
steep power-law one, which is consistent with our model. 

What should be investigated for a further understanding of flux transport may be a 
statistical investigation of cancellation on the photosphere. Another important question is 
whether there is a stable solution of our model with four magnetic processes or not. We 
will investigate these questions in the series of this paper. 
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2009 November 11 In 2:21 UT 




-40 -20 20 40 
Solar-X (arcsec) 

Fig. 1. — An example of Na I D 1 magnetograms after the pre-process explained in the 
text. This magnetogram is taken at 2:21UT on 2009 November 11, which is the middle of 
observational period. Some magnetic networks are contained in the field of view. 
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SOT CP vs MDI Magnetic Field 
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Fig. 2. — A scatter plot of CP obtained by Hinode/SOT and magnetic flux obtained by 
SoHO/MDI. Asterisks show averaged SOT CP signal corresponding to one pixel of MDI. 
The bars in horizontal direction indicate minimum and maximum value of SOT CP signal 
in each MDI pixel. Solid line shows a result of a linear fitting, whose slope is 9067.38 G/DN. 
The fitting range is from 30 G to 100 G in absolute magnetic field strength. 
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Fig. 3. — Schematic pictures of four magnetic processes between two patches, (a) Merging. 
Two patches converge and coalesce into one patch, (b) Splitting. One patch is divided into 
two patches, (c) Emergence. Two opposite polarities with equal amount of flux content 
appear, (d) Cancellation. Two opposite polarities converge with each other. One of them 
disappears. The other loses its flux content or disappear. 
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Fig. 4.— Schematic picture of our patch tracking and magnetic process detection. 1. 
Detection of patches with clumping algorithm. 2. Track of patches by examining overlaps. 
3. Detection of merging and splitting by examining overlaps. 4. Detection of cancellation 
and emergence by making pairs of flux change events. 
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(c) S th =81 (pixel) (d) S th =144 (pixel) 
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Fig. 5. — Two-leveled magnetograms taken at 2:21 on 2009 November 11 with different size 
thresholds. Size thresholds, S th , are set as (a) 9 pixels, (b) 36 pixels, (c) 81 pixels, and (d) 
144 pixels. 
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Fig. 6.— Examples of detected magnetic processes, namely (a) merging, (b) splitting, 
(c) emergence, and (d) cancellation. Background shows magnetic flux density obtained by 
Hinode/NFl. Red (blue) contours indicate positive (negative) patches detected with our 
threshold. The field of view is 14.4" x 14.4" for all the images. 
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Frequency Distribution of Flux Content 
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Fig. 7. — Frequency distribution of magnetic flux content. Red (blue) dashed histogram 
denotes that of positive (negative) polarity. Solid histogram denotes sum of them. Black 
dashed line indicates the fitting result of it. The obtained power-law index is -1.78 with a 
fitting; ranee between 10 175 Mx and 10 19 Mx. The vertical green line denotes the detection 
limit, 0th= 10 17 - 5 Mx. 
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Fig. 8. — Apparent probability distribution of merging. It is made from merging where the 
parent patches have flux content more than (p t h- Red/blue/black solid lines indicate obser- 
vational results for positive/negative/both patches. Red/blue/black dashed lines indicate 
fitting results with a range of 10 17 ' 5 - 10 19 Mx. The power-law indexes of the fitting lines are 
0.28, 0.26, and 0.28. Horizontal and vertical dashed line indicates to timescale of 33 minutes 
and the detection limit, respectively. 
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Fig. 9. — Apparent probability distribution of splitting. Red/blue/black solid lines indicate 
observational results for posi- tive/negative/both patches. Blue and black dashed lines in- 
dicate analytical curves with ko = 5.0xl0~ 4 s _1 and 1.0 xlO -3 s -1 . Horizontal dashed line 
indicates a timescale of 33 minutes. 
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Fig. 10. — Model of relationship among flux distributions of flux content, cancellation, and 
emergence with 1.5 < 7 < 2. Solid arrows indicate mathematical relations, while dashed 
arrows indicate relationships with some assumptions. 
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Table 1: Total numbers and flux amounts of magnetic patches and magnetic processes. The 
numbers in the parentheses in the first line are averaged flux densities of each polarity. Those 
in the parentheses in the second and third lines are averaged rates of flux density change 
by each process. These values are averaged over the whole field of view and observational 
period. 



